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We investigate long-range forces between atoms with anisotropic electric polarizability interacting 
with dielectrics having anisotropic permittivity in the weak-coupling approximation. Unstable con- 
figurations in which the force between the objects is repulsive are constructed. Such configurations 
exist for three anisotropic atoms as well as for an anisotropic atom above a dielectric plate with a 
hole whose permittivity is anisotropic. Apart from the absolute magnitude of the force, the depen- 
dence on the configuration is qualitatively the same as for metallic objects for which the anisotropy 
is a purely geometric effect. In the weak limit closed analytic expressions for rather complicated 
configurations are obtained. The non-monotonic dependence of the interaction energy on separation 
is related to the fact that the electromagnetic Green's dyadic is not positive definite. The analysis 
in the weak limit is found to also semi-quantitatively explain the dependence of Casimir forces on 
the orientation of anisotropic dielectrics observed experimentally. Contrary to the scalar case, ir- 
reducible electromagnetic three-body energies can change sign. We trace this to the fact that the 
electromagnetic Green's dyadic is not positive definite. 

PACS numbers: 31.30.Jv, 32.10.Dk, 33.15.Dj, 34.20.Gj, 77.22.-d, 81.05.Xj 



I. INTRODUCTION 

Fifty-seven years after the postulation of attractive long-range van der Waals forces between neutral atoms [| 
London described them as a purely quantum mechanical effect. Later Axilrod and Teller and Muto 

independently found that the three-body contribution of three atoms to the van der Waals-London interaction energy 
can change sign and contribute repulsively depending on their configuration. Repulsive long-range forces have so far 
been established in at least two situation£0: for two dielectric objects with a suitable dielectric between them 0, Q , 
and for a dielectric object and a magnetic material IQ., JdJ. The latter repulsive effect has not yet been observed. 
Although it was known that the interaction between anisotropic atoms could have repulsive components [H, , it 
came as a surprise that the force on a vertically oriented metallic needle above a hole in a metallic plate can be 
repulsive as well [l^ . The explanation of this effect relied on symmetry arguments and numerical calculations. An 
analytic solution, in the unretarded regime, for a polarizable particle above a plate with a hole was obtained in [isj 
with similar repulsive effects. By summing two-body contributions in the multiple scattering formalism the effect 
could also be verified in the retarded regime The origin of such repulsive effects nevertheless does not seem to 
have been sufficiently well understood to make reliable predictions for more general situations. For instance, although 
very interesting repulsive configurations were found in (TtI . [18| one of the conclusion was that "there is no repulsion 
possible in the weak coupling regime" . For anisotropic atoms and dielectrics this is not the case. 

The interaction between anisotropic atoms has been studied extensively (see [IS] and references therein). We in 
the following identify configurations with repulsive long-range forces between atoms and dielectric materials with 
anisotropic polarizabilities and permittivities. For weak polarizability and susceptibility we obtain exact analytic 
expressions for a wide range of geometries. The force on an atom sufficiently close to a dielectric plate with a hole 
is found to be repulsive for certain relative orientations of the polarizability and permeability of the atom and the 
dielectric material. It is always attractive at large separations. The cause for Casimir repulsion is traced to anisotropies 
in the polarizability of the objects and to the fact that the electromagnetic Green's dyadic is not positive definite. 
For ideal metals the anisotropy is geometrical and exact results for dilute dielectrics should qualitatively extend to 
metals. 

We begin by examining long-range forces between anisotropic atoms and dielectrics. As for scalar Casimir energies 
[20l - |23j we can show that the two-body contribution to the electromagnetic Casimir interaction is always negative, 
independent of the relative orientation of polarizability and permittivity tensors. However, contrary to the scalar 
case, the electromagnetic Casimir interaction energy in general is not a monotonic function of separation and some 
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^ We here do not include the repulsive pressure on a sphere found by Boyer jTJ because it does not involve two objects. 
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force-components can change sign. Remarkably, torque-free points at which the interaction energy does not depend 
on the orientation of the atom appear to exist for some configurations of dielectrics. We then examine the three-body 
correction to the interaction of isotropic atoms [131 and find that it never dominates two-body contributions to 
the Casimir energy. Contrary to the scalar case, irreducible three-body contributions to the electromagnetic Casimir 
energy can change sign, because the electromagnetic Green's dyadic is not positive definite. 



II. ANISOTROPIC POLARIZABLE ATOMS 



Long-range interactions of polarizable atoms are described by the van der Waals-London forces [1, Q in the un- 
retarded regime and by Casimir-Polder forces [25| in the retarded regime. In the multiple scattering formalism the 
interaction energy for two polarizable atoms is 
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an expression already obtained for atoms in [2^. The free Green's dyadic Tq in Eq. ([T]) is 
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with u{x) = 1 + X + x'^, and v{x) = 3 -f 3a; -I- x'^. Neglecting quadruple and higher moments, the scattering matrix Ti 
for the i-th atom with atomic dipole polarizability Q!i(iC) is 



T,(x,x'; iC) = 47ra,(iC) (5(3)(x - xJ^J^'^x' - x,), 



(3) 



where x^ specifies the position of the atom. The (5-functions in Eq. ^ permit a trivial evaluation of the spatial 
integrals of the trace in Eq. ([1]). For separations r^^ satisfying rfj |ai(0)aj(0)| the logarithm in Eq. ([IJ may be 
expandecQ. The weak approximation consists of retaining only the leading term of this expansion and we have 
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which is the interaction obtained in [1^ |13|. We used the property that polarizability is a symmetric tensor, of the 
form (Xi = ^nCt^e^e"^ , where the e", n = 1,2,3, are the orthogonal principal axes satisfying e™ • e" = S"^" and 
e™ X e" = e™"'e'. For stable atoms the corresponding principal polarizabilities a" necessarily are non- negative. A 
simple model for the frequency dependence of the atomic polarizability is 



(5) 



where Wi is the excitation energy of the (two-level) atom and ai{0) is its static polarizability. One of the ways to 
calculate atomic polarizability is by Dalgarno's method [271] . 

The exponential dependence on the separation distance in Eq. ^ implies that the frequency dependence of the 
polarizability ai{iC,) in Eq. ([5|) is negligible for < c/uji <C r. In this asymptotic retarded (Casimir-Polder) regime 
the polarizabilities can be approximated by their static values and the ^-integration in Eq. ^ performed to yield 

[a El, 



^^12 ("i,"2;r) 



87rr7 



13Tr{ai(0) •Q:2(0)}-56{f •ai(0) •a2(0) •f}-|-63{f -0:1(0) -fllf -0:2(0) -f} . (6) 



For atoms with isotropic polarizabilities, ai = ail and 02 = 0:21, Eq- © reproduces the Casimir-Polder interaction 
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^ Approximating atoms as point-like objects in Eq. ((Sjl is not justified for rf- < \ai{0)cs.j(0)\. 
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FIG. 1. Two atoms with anisotropic polarizabilities. 



In the unretarded (London) regime {a^ < r <^ c/uii), the frequency dependence in Eq. ([2]) may be neglected and 
the free dyadic in Eq. ^ approximated by the static dipole-dipole interaction ro(r;0) to yield [l2l. [Tsj. 
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For atoms with isotropic polarizabilities this reproduces London's expression for the van der Waals interaction 
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which is inversely proportional to the sixth power in the separation r. To evaluate the coefficient, the frequency- 
dependence of the polarizabilities has to be known or modeled. For the simple model in Eq. ([5]) (letting wi ^ uj2 = t^o) 
one has 
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With the increased computational power it nowadays is only slightly more complicated to perform the complete 
frequency integral in Eq. (j4]) numerically. 

To illustrate the orientation dependence of the Casimir-Polder interaction, consider the special case where one of 
the atoms is isotropic and the other a linear molecule that essentially can only be polarized along its axis with a 
polarizability given byQ;i=aiee. With e • f = cos 9, Eq. (jH]) then becomes 
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The minimum energy configuration is at = or = tt, corresponding to the alignment of the axis e of the molecule 
with r. Note that the interaction energy of Eq. (1111) is negative for any orientation. 

Consider next two atoms with general static anisotropic polarizability tensors = J2n'^7^?^2 ■ (^^e Fig. [1]) 
Inserting these in Eq. ([SJ, the interaction energy in the retarded regime is of the form 



i?pP(ai,a2;r) = --i^ ^ aj" (0)a^(0)CiT, 
onr' ^ — ' 
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with 
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The potential energy surface of the interaction given in Eq. (|12p is rich, with (local) extrema of individual C{^" at 
C™^ = 20,13,5/9, and 0. These extremums are summarized in Table H] The C"2"'s are non-negative and vanish 
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±1 
±1 
±1 





±1 
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-20 



-13 



Minima 
Saddle 
Saddle 

Maxima 




TABLE I. Configurations with extremal energies for two linearly polarizable molecules. In the rows corresponding to C'Yi" = 20, 
and 5/9, the signs in the first two columns should concur. 



only if e^" • e" = and (e^" • f)(e2 • f) = 0. Since the eigenvalues of the polarization tensor are non- negative, the 
two-body interaction between two anisotropic atoms therefore is always negative. The maximum value of C{^" = 20 
corresponds to the energetically most favorable orientation, e™ = ±r, = ±r, where both atoms are aligned with 
r. The potential energy surface has a saddle point when (e™ • f)(e2 • r) = and e™ • e" = ±1 that corresponds to 
C{^" = 13. Interesting extremum for (e™ • f) = ±2/3 = (ej ■ r) and e™ • e" = ±1 corresponds to = 5/9. Note 
that (e™ • f) = ±2/3 represents an angle of about 48.2° between the polarizabilities and f. The interaction energy in 
Eq. (jl2p gives rise to a non-central force between anisotropic atoms, 
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The second term in Eq. (fT4|) is a torque that vanishes for the extremal configurations of Table HI 



III. CASIMIR REPULSION 



Since the Casimir-Polder energy can vanish for particular orientations of anisotropic atoms and is always negative, 
it is clear from the foregoing discussion that it in general is not a monotonic function of the distance between the 
atoms and that for fixed orientation of the atoms components of the force between them can be repulsive. To see 
this more explicitly consider two linear atoms with polarizabilities that are orthogonal to each other; ai — aiaa, 
a.2 = a2hh, with a • h = 0. Let a • f = cos 6*, h • f = sinO, = ± h^. (See insert in Fig. [H) Eq. ^ gives the 
interaction energy for this configuration as 
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For fixed horizontal separation a, the dimensionless interaction energy of Eq. (jl5p is shown in Fig. [51 It vanishes for 
h = and h oo and is negative otherwise. The minimum at h — h/a — •\/2/3 ^ 0.47 implies that the vertical 
component of the force on atom 2 along h changes sign when 3h = \f2a. Regimes with repulsive components of the 
non-central force between anisotropic atoms exist for all values of a • h. For |a • h| sufficiently close to unity, the 
repulsive regime splits into two or more disjoint regions in h. In particular, for |a • h| = 1 the two repulsive regimes 
are disjoint, symmetrical, and above and below the plate. 

In the previous example, the horizontal component of the force between the two atoms was always attractive. 
However, this component of the force on atom 2 may be compensated by placing a third atom with polarization 
tensor 0:3 = ±aiaa in the symmetric point with the three atoms forming an isosceles triangle. Of course, real atoms 
would not maintain these orientations and positions unless they are part of a material. These considerations bring 
us to the analog of the geometry of a metallic needle above a metallic plate [13] • We therefore consider an atom with 
anisotropic polarizability centered above a dilute dielectric material with a hole as sketched in the inserts of Fig. [S] In 
the weak limit the scattering matrix of dilute dielectric objects can be approximated by their dielectric permittivity 
with Ti(x, y) ^ V^(x)5^^^ (x — y) and Vi(x) = [e(x) — 1] ^ 47rQ:i7T,i(x), where ni(x) is the number density of the atoms 
the dielectric is composed of and on is their polarizability. In this weak limit Eq. ([T]) in the retarded regime takes the 
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FIG. 2. Dimensionless energy, with Eo = 63ai (0)Q2(0)/87r a^, for two atoms with orthogonal polarizabilities in Eq. (|15|l . 
sketched in the insert, plotted with respect to h/a. 
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Eq. (|T6|) generalizes the corresponding expression for isotropic dielectric functions in ^] to the anisotropic case. 

Using Eq. (fTB]) we now can study the analog of the configuration proposed in [3] for dilute anisotropic dielectrics 
and atoms. In a Cartesian coordinate system with axes in the x, y, and z directions, we consider a dielectric plate 
in the xy-plane with a hole of radius a centered at the origin and an anisotropic (linear) atom at the height h on 
the z-axis. The scattering potential of the anisotropic atom is modeled by, Vi(x) = ATTa.5{x)5{y)5{z — h) with the 
anisotropic static polarizability a = a e e and e i — cos 9. The scattering potential of the dielectric plate of thickness 
d with a hole of radius a is described by 



V2(x) = (£ - 1) e{p - a)[e{z + d)- e{z)] ^\e{p- a)5{z), 



(17) 



with — x"^ + . To simplify the calculation we consider a thin plate and define A ^ (e — 'l)d ~ Airaa, with a = nd 
the planar density of atoms. We furthermore assume that the dielectric material is polarizable in the plane of the 
dielectric only and therefore demand that A • z = 0. If A is isotropic in the xy plane we can insert A = A (xx -t- yy) 
in Eq. ([TS]) to obtain 
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for the interaction energy in the retarded regime. Fig. [3] shows the dependence of Eq. ([T5| on the dimensionless 
height h ^ h/a for different orientations 9 of the anisotropic atom. At the intersection of all the curves in Fig. [3] 
the interaction energy does not depend on the orientation 6 of the atom and the anisotropic atom is torsion free. 
This is the case when the coefficient of cos 20 in Eq. ((T5)) vanishes, that is for 4ft,'' + 83h^ — 26 = 0, giving h ^ 0.56. 
The torsion-free point coincides with the highest (unstable) equilibrium point, attained for 9 — 0. For heights above 
the torsion-free point the energy is minimized for vertical orientation of the polarizability and below it the energy is 
minimized for horizontal orientations. The transition from vertical to horizontal orientation is sudden and happens 
at the torsion-free point without expense in energy. The net force on the atom in z-dircction is 
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and is shown in Fig. [3] as a function of h for several orientations 9 of the atom. The orientation-independent value for 
the vertical force occurs when Ah"^ + 113ft-^ — 80 = and corresponds to h ^ 0.83. The force on the atom given by 
Eq. (|19p is repulsive when 
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FIG. 3. The dimensionless energy in Eq. (|18|) with Eq = aA/3207ra^ (left) and the dimensionless force in Eq. H19[) with 
-Fo = aX/GAna^ (right) as functions of h/a for an anisotropic atom at a height h above the center of a dielectric plate with a 
hole of radius a as sketched in the inserts. The curves correspond to different orientation 9 of the atomic polarizability. 



which has real solutions for 

0<9<lcos-'(J^^ and n - ^cos-' (J^^ < 9 < tt. (21) 

The critical value for the orientation angle, 9c = 0.5cos~^(l/10) ~ 42.1°, is the angle beyond which no repulsive 
regime exists. 

To emphasize how weakly interacting configurations could lead to a semi-quantitative understanding of their strongly 
interacting counterparts, we next consider the configuration in [295, where the Casimir force between gold and the 
anisotropic cuprate superconductor Bi2Sr2CaCu208+(5 (BSCCO) was found to decrease in magnitude by as much as 
10-20% when the BSCCO planes were oriented parallel to the gold plate rather than perpendicular to it. Let us 
replace the gold plate by a semi-infinite dielectric slab with isotropic permittivity described by Vi(x) — [ei — l]d{~z) 
with (ei — 1) = (ei — 1)1, and replace the BSCCO by a semi-infinite dielectric slab with V2(x) — [e2 — i]0{z — a) 
and (£2 — 1) = (£2 — 1)± (xx + yy) + (£2 — l)||zz when the BSCCO-planes are parallel to the xy-plane (perpendicular 
cleave in [29'|) and (£2 — 1) = (£2 — l)||xx -I- (£2 — 1)± (yy -I- zz) if they are oriented parallel to the xz-plane (parallel 
cleave in [29,]). Using Eq. one obtains for the ratio of the Casimir force in the two orientations 

■^11 clcavo ^ ^ 2 [(£2 - 1)± - (£2 - 1)||] 
cleave [l4(£2-l)±+9(£2-l)||]' 

If we assume (£2 — 1)|| <C (£2 — 1)_l, which is a good approximation for BSCCO, the magnitude of the Casimir force 
in the retarded regime changes by 1/7 ~ 14%, in semi-quantitative agreement with [29j . 



IV. THREE-ATOM CASIMIR ENERGIES 

Irreducible many-body Casimir energies were proven to remain finite in pol [23j when some, but not all, objects 
overlap. For a massless scalar quantum field with local potential interactions, the sign of the A^-body Casimir energy 
was found to be simply (—1)^+^. These theorems were verified in all examples studied in [20l - |23| where closed 
expressions for irreducible many-body Casimir energies were derived in the framework of the multiple scattering 
expansion. The proof of the sign of irreducible many-body contributions in the scalar case relied on the positivity of 
the free scalar Green's function. The electromagnetic free Green's dyadic of Eq. 1^ is not a positive definite operator. 
The eigenvalues of this matrix, corresponding to eigenvectors parallel to r and orthogonal to it, are proportional to 
(^(a;) — v{x),u{x),u{x)), where u(x) and v{x) were defined after Eq. Here we investigate consequences of this for 
the sign of the irreducible three-body contribution to the Casimir energy between three isotropic atoms 0, . 

Inserting the free Green's dyadic of Eq. ([2]) into the expression for the irreducible three-body contribution to the 
Casimir energy in and using scattering matrices for isotropic atoms given by Eq. with ai = ai{i()l one 
obtains 

1 1 r"^ 

^1^3 - - 3 3 3 / c;Ce-«('--+'--+'--)ai(zC)a2(*C)a3«)Qi23(ri2,r23), (23) 

^ ^12^23' 31 Jo 
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FIG. 4. Angular dependence of three-body contributions to the long-range forces in the retarded (CP) and unretarded (Lon) 
regime for three atoms forming an isosceles triangle. 



where the displacements Vij are the sides of the triangle formed by the three atoms and satisfy ri2 + + r^i = 0. 
The function Q123 in Eq. is 



(24) 



Ql23(ri2, r23) = 3M12U23M3I - W12U23W3I - U12W23U3I - t'l2U23U31 + Wi2t;23W31 COS^ 63 
+U12U23W3I COS^ 9i + ^;12^'23U31 COS^ 62 + Wl2W23t'31 COS 0i COS 6*2 COS 03, 



with Uij — u{\C\rij), Vij — w(|Ckii)i the Si's (z — 1,2,3) being the internal angles of the triangle. In the unretarded 
limit Eq. (|23p simplifies to the expression obtained in 15], 



Lon _ (1 +3C0S6>1 COS6I2 cos6>3) 3 
''^123 — „3 „3 „3 ^ 
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Using the simple model in Eq. (I5|) for three identical isotropic atoms with uji = uj2 = W3 ~ luq, the integral in Eq. ((25 
is readily performed to yield ^30.] . 
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where the angular dependence is given by |24| 

(^123(^1 ^^■2) 



9 sin^ ^3(1 + 3 cos 9i cos 6*2 cos 63) 
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(27) 



with the three sides and angles of the triangle related by the law of sines. It was noted in [5] that the three-body 
contribution in the unretarded regime is negative when the atoms form an acute triangle and positive when it is 
(very) obtuse. For isosceles triangles this behavior is seen in Fig. 01 Refs. (30l - [3^ attempted to explain the crystal 
structures of rare gases by three-body contributions to the total energy. Although this turned out not to be possible, 
three-body effects typically contribute between 10-20% to the total energy. 

The C-integral in Eq. (j23p can be performed analytically in the retarded limit but the result is not very illuminating 
and too lengthy to be presented here. In the special case when the atoms are arranged to form an isosceles triangle 
with side lengths ri2 = a, r23 = r^i — r and 9i = 62 ~ 0, (see insert in Fig. U) the retarded interaction is [23] 
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7 -I- 49a; -I- eila;^ + 1533x^ + 868a:'' - 1372a;^ - 1672a;^ - 672a;^ - 96x^ 



(28) 



(29) 



For an equilateral triangle it yields Ey^^{a,TT/3) — ai (0)Q;2(0)Q;3(0)1264/2437ra^° in agreement with [2^]. As in the 
unretarded limit, the angular dependence in the retarded limit also leads to configurations in which the three-body 
force contributes repulsively. Fig.|4]in fact shows that the irreducible three-body contribution to the potential is less 
attractive in the retarded regime. However, for isotropic atoms, the three-body contribution never dominates over 
two-body contributions in the regime a < where a point-like description of the atoms is justified. Evidently, the 
analysis of atomic many-body interaction is only a little more involved for anisotropic atoms and molecules. 
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V. CONCLUSION 



We have shown that the force between atoms with anisotropic polarizabiHties and dilute dielectrics with anisotropic 
permeabilities can have repulsive components. The weakly interacting configurations we considered give a semi- 
quantitative understanding of the analogous strongly interacting case. That the two-body Casimir energy is not always 
monotonic in the separation is associated to the tensorial structure of the polarizabilities and the electromagnetic 
Green's dyadic. The analysis of an anisotropic atom above a dilute dielectric plate with a hole provides considerable 
insight into the analogous configuration involving perfect metals considered in p^ . Closed analytic solution for this 
configuration reveals torque- free points at which the interaction energy is independent of the orientation of the atomic 
polarizability. Although a single hole in a dielectric plate does not lead to stable points, multiple holes are expected 
to introduce stable points in the potential energy surface of an oriented atom. 

Unlike the scalar case, for which the irreducible three-body contribution to Casimir energy is always positive (20l - [23| . 
we find the electromagnetic three-body Casimir energy can change sign. The three-body contributions to the energy 
are found to never dominate the two-body contributions for isotropic atoms. 
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